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ABSTRACT 


A comparison  of  two  air-sea  interaction  measurement  systems  was 
run  at  Gable  Island,  N.  S.,  in  September  and  October,  1975.  The  BIO 
sensors  included  thrust  and  sonic  anemometers  for  wind  turbulence 
measurement,  a microbead  thermistor,  and  a Lyman-alpha  humidity  sensor. 
TheJJBC  system  consisted  of  a modified  Gill  propeller -vane  anemometer, 
^'"thermistors,  and  Brady  and  aluminum  oxide  type  humidity  sensors.  Both 
data  sets  were  recorded  in  fm  analogue  form  and  digitized  and  analyzed 
at  BIO.  The  UBC  system  also  made  digital  recordings,  which  will  be 

reported  elsewhere. 
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?Both  systems  operated  together  to  produce  generally  comparable 
results.  A certain  number  of  differences  and  problems  are  discussed. 


SOMMAIRE 


En  September  et  octobre  1975 » A l’ile  de  Sable  (N.-E.),  on  a 
proc£d£  a la  comparaison  de  deux  systAmes  de  mesure  de  1 ' interaction 
atmosphere-mer . Les  dAtecteurs  de  l'Institut  ocAanographique  de  Bedford 
comprenaient  des  an^mom&tres  soniques  et  de  poussAe  destines  a mesurer 
la  turbulence  des  vents,  un  micro-thermistor  et  un  dAtecteur  d' humidity 
Lyman-A 1 phn.  Le  systeme  de  1 'University  de  Colombie-britannique  se 
composaLt  d'un  an^mom&tre  Gill  modifie,  a h^lice  et  girouette,  de 
thermistors,  et  de  d^tecteurs  d'humidite  de  types  Brady  et  A oxyde 
d'aluminium.  Les  deux  series  de  donnAes  ont  enregistrAes  de  faqon 
analogue,  en  FM,  et  analyses  A l'Institut  oc^anographique  de  Bedford. 

Le  systAme  de  1 'University  de  Colombie-britannique  a aussi  fait  des 
enregistrements  num^riques  dont  le  rapport  figurera  ailleurs. 

Les  deux  systAmes  ont  fonctionnA  ensemble  pour  donner  des 
resultats  comparables  dans  1' ensemble.  II  est  discutA  dans  le  texte 
d'un  certain  nombre  de  differences  et  de  problAmes. 
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SYMBOLS 

* 

B0’WB3 

C10=-“lu3/U102 

Coefficients  of  polynomial  relating  thrust 
anemometer  signals  to  wind  tilt  in  degrees. 
Sea  surface  drag  coefficient. 

Heat  flux  coefficient 

WW« 

Evaporation  coefficient. 

f 

Frequency 

GMT 

Greenwich  Mean  Time,  4 hours  ahead  of  local 

i 

Atlantic  Standard  Time. 

Q 

Mean  water  vapour  density 

ij 

i 

q 

rl3=  uiU3/o1°3 

Fluctuation  in  water  vapour  density. 
Correlation  coefficient  Between  downwind  and 
vertical  wind  components. 

T 

Mean  air  temperature 

i 

t 

Air  temperature  fluctuation 

T 

s 

Sea  surface  ("tucket")  temperature. 

AT=T-T 

Temperature  difference  between  sea  surface  and 

4 

10  m height. 

U 

Mean  wind  velocity. 

• 

uio 

Mean  wind  velocity  at  a height  of  10  m above 

the  surface 

U. 

1 

Wind  velocity  component. 

3. 

u#=(-u1u3)s 

Friction  velocity. 

VH 

Horizontal  velocity,  in  anemometer  axes 

- 

W 

Vertical  velocity,  in  anemometer  axes 

i 

W 

vS  >a 


u.  ,1)  , 

t q 


■i=l,2,3(as  subscript) 


overbar 


IV 


Standard  deviation  of  velocity  component. 


Standard  deviation  of  temperature,  humidity. 


Spectrum  (if  i=j;,  or  cospectrum  (if  i^j ) 


Downwind,  lateral,  and  vertical  components. 

(The  first  axis  is  aligned  with  the  mean  wind.) 


Average  over  a data  run,  typically  of  about 
40  min  duration. 


1.  INTRODUCTION 

Air-sea  interaction  measuring  systems  have  been  designed  and  con- 
structed at  the  Bedford  Institute  of  Oceanography  (BIO)  and  at  the  Institute 
of  Oceanography,  University  of  British  Columbia  (IOUBC)  to  measure  wind 
velocities,  humidity  and  temperature  for  determinations  of  momentum,  moisture 
and  heat  fluxes  between  the  ocean  and  atmosphere.  Although  the  two  systems 
are  intended  to  measure  the  same  physical  quantities,  they  use  different 
sensors  and  data  acquisition  systems.  In  addition,  the  analogue  tape 
recorder  in  the  BIO  system  has  relatively  high  power  consumption  whereas 
the  IOUBC  system  has  very  modest  power  requirements.  Further,  the  IOUBC 
system  records  digitally  in  intermittent  bursts  to  greatly  reduce  the 
quantity  of  tape  required  with  the  result  that  it  can  be  self-contained, 
while  for  remote  sites  the  BIO  data  must  be  recorded  using  telemetry. 

An  intercomparison  experiment  was  conducted  with  the  two  systems  in 
order  to  be  able  to  compare  future  results,  especially  those  to  be  obtained 
from  a stable  tower  moored  off  Halifax  Harbour.  In  the  present  experi- 
ment, data  from  both  the  BIO  and  IOUBC  sensors  were  recorded  on  the  BIO 
data  recording  system.  The  IOUBC  system  records  will  be  compared  to  these 
results  when  the  IOUBC  analysis  system  development  is  completed. 

The  intercomparison  study  was  made  on  Sable  Island,  Nova  Scotia, 
from  mid-September  to  the  end  of  October  1975.  Sable  Island  was  chosen 
because  it  is  a relatively  convenient  site  for  air-sea  interaction  studies 
and  because  earlier  work  had  been  done  there  with  the  BIO  system  (Smith  and 
Banke,  1975).  This  earlier  work  was  done  on  the  west  spit  of  the  island, 
but,  because  reaching  that  site  from  the  support  facilities  of  the  meteoro- 
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logical  station  was  time  consuming  and  sometimes  impossible  due  to  beach 
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Figure  1.  The  instrument  tower  at  low  tide.  The  high  water  line  can  be 
seen  half  way  up  the  beach.  Inset  - the  sensors:  (A)  UBC 

Gill  anemometer;  (B)  shield  containing  UBC  temperature  and 
humidity  sensors;  (C)  BIO  Mark  7 thrust  anemometer; 

(D)  mounting  bracket  for  BIO  Lyman-alpha  humidiometer ; 

(E)  BIO  sonic  anemometer  with  thermistor  attached. 


conditions,  a more  convenient  site  on  the  south  side  of  the  island  nearer 
the  meteorological  station  was  chosen  for  this  work. 


2.  INSTRUMENTATION 

The  various  sensors  were  mounted  10  to  13  m above  the  sea  surface 
(depending  on  the  tide)  on  a 9 m high  mast  located  about  30  m from  the 
water's  edge  (Fig.  l).  The  IOUBC  digital  recording  system  was  mounted  at 
the  base  of  the  mast  and  the  BIO  recording  system  was  housed  in  a hut  about 
50  m inland  from  the  mast. 
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The  BIO  sensors  consisted  of  a Kaijo-Denki  PAT  311-1  sonic  anemom- 
eter, two  thrust  anemometers,  a Victory  Engineering  Corp.  microbead 
thermistor,  and  an  Electromagnetic  Research  Co.  Model  BLR  Lyman-alpha 
humidity  sensor. 

The  Kaijo-Denki  three-component  sonic  anemometer  has  been  widely 
accepted  for  measurements  of  atmospheric  turbulence  (Japan-U.S.  Joint  Study 
Group,  1971)  and  served  in  this  experiment  as  a basis  for  comparison  of  the 
other  anemometers.  This  anemometer  has  one  failing:  when  the  probes 

become  wet  from  rain,  the  sound  pulses  are  not  always  detected  properly  and 
the  signals  contain  large  spikes.  In  addition  it  requires  frequent 
adjustment  and  consumes  about  100  W,  and  so  will  not  be  used  on  the 
stable  tower. 

The  two  thrust  anemometers  measure  three  components  of  wind  thrust 
on  a perforated  ping-pong  ball  using  spring  linkages  and  differential 
transformers.  The  Mark  6 design  (Smith,  1969)  has  been  in  use  for  some 
time.  its  spring  linkage  suffers  from  thermal  drift,  and  from  cross- 
eoupJ ing  of  the  three  axes  which  introduces  errors  up  to  3°  in  indicated  wind 
direction.  The  anemometer  used  was  the  fourth  one  built  to  this  design  and 
will  be  referred  to  as  Thrust  6.U.  Thermal  drift  can  be  essentially  elim- 
inated by  recording  zero-wind  signals  immediately  before  and  after  each 
data  run  and  by  the  removal  of  linear  trends  from  the  data.  Some  of  the 
cross-coupling  errors  can  be  compensated  for  in  the  data  processing,  but 
variations  in  cross  coupling  with  wind  thrust  and  direction  cannot  be 
eliminated  entirely.  The  Mark  7 design,  referred  to  as  Thrust  7.1,  has 
a simplified  spring  linkage  (Fig.  2)  with  reduced  cross-coupling  intro- 
ducing a smaller  error,  about  1°,  in  the  indicated  wind  direction.  It 
lacks  the  counterbalancing  against  horizontal  acceleration  of  the  Mark  6 


and  so  is  more  prone  to  errors  caused  by  movements  or  vibrations  of  its 
supporting  mast.  Because  of  this  problem,  which  was  very  apparent  in 


the  Sable  Island  measurements,  the  development  of  the  Mark  7 anemometer 
has  been  discontinued.  It  remains  a usable  instrument  if  mounted  on  a 
very  solid  tower. 

The  BIO  microbead  thermistor  (Thorpe,  1973)  was  exposed  to  the  air 
stream  in  a small  (25  mm  diameter)  radiation  shield  mounted  near  the  sonic 
anemometer  probe. 

The  Lyman-alpha  humidiometer  was  operated  with  an  absorbing  path 
length  of  0.5  cm  and  a source  current  of  50  yA.  Calibration  was  achieved  by 
placing  it  in  an  enclosed  box  with  a psychrometer  and  a Cambridge  Systems 
Model  137  dew  point  hygrometer.  A number  of  calibration  points  over  a small 
range  of  absolute  humidity  was  obtained  in  this  manner  and  the  ambient 
humidity  reading  was  checked  with  a psychrometer  on  the  tower  before  each 
data  run.  Some  drift  did  occur.  Because  the  humidiometer  was  exposed  to 
the  air  flow  with  no  weather  shield  and  was  vulnerable  to  moisture  damage, 
it  was  mounted  on  the  tower  only  during  some  data  runs. 

The  lUO  data  were  recorded  on  IRIG  (intermediate  Range  Instrumenta- 
tion Group ) frequency-multiplexed  channels  1 to  12  using  a Crown  Model  Cl 
821’  tape  recorder  running  at  V-i  inch/s.  A 1^.5  kHz  reference  tone  was  added 
to  the  signal  to  allow  compensation  during  playback  for  tape  speed  variation. 

The  I0UBC  system  measured  wind  velocities  with  a modified  Gill 
anemometer  (R.M.  Young  & Co.,  Traverse  City,  Michigan,  U.S.A.)  which  has 


the  standard  propeller  vane  unit  and  a lengthened  vertical  shaft.  A second 
propeller  was  added  with  its  axis  pointing  downward  at  an  angle  of  about 
30°  from  the  vertical  instead  of  pointing  vertically  upward  as  is  more 
usual.  The  non-vertical  axis  avoids  the  stalling  and  non-linear  response 
common  in  the  usual  configuration  for  vertical  wind  velocities  less  than 
H m/s;  The  tilted  propeller  ensures  that  it  always  operates  above  the 
threshold  and  non-linear  regions  for  wind  speeds  above  3 m/s.  The 
second  propeller  thus  measures  a combination  of  the  horizontal  component  of 
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the  wind  velocity,  V , and  the  vertical  component,  W,  but  it  does  not 

n 

follow  a cosine  law  exactly  and  a more  complicated  calculation  is 
required  to  extract  W from  the  two  propellers. 

The  IOUBC  system  also  has  a microbead  thermistor  (Victory  Engin- 
eering Corp. ) for  temperature  measurements,  a glass  rod  thermistor  to  meas- 
ure sea  temperature  which  in  this  experiment  was  used  to  check  the  micro- 
bead thermistor,  a Brady  array  (Thunder  Scientific  Inc.  ) for  relative 
humidity  measurements  and  an  aluminum  oxide  sensor  (Panametrics  Inc. ) for 
absolute  humidity  measurements.  The  temperature  and  humidity  sensors  were 
mounted  in  a protective  shield  of  the  size  and  shape  of  a sombrero  (inset 
in  Fig.  1).  The  IOUBC  sensors,  the  method  to  obtain  u^  from  the  Gill 
anemometer,  the  IOUBC  dissipation  method  and  eddy-correlation  method 
recording  packages  will  be  described  elsewhere  (Pond  and  Large,  1977). 

3.  LATA  PROCESSING 

The  analogue  data  tapes  were  digitized  at  the  BIO  Computing  Centre 
using  Program  A TO  D (Smith,  197M-  This  program  linearized  the  thrust 
anemometer  and  BIO  thermistor  output,  converted  the  sonic  anemometer  data 
to  orthogonal  coordinates , and  used  range  and  offset  values  to  convert  data 
from  linear  sensors  to  physical  units.  In  a second  pass  through  the  data 
the  program  removed  means  and  linear  trends,  rotated  coordinates  from  the 
anemometers  into  coordinates  aligned  with  the  indicated  mean  wind  direction 
so  that  the  vertical  and  cross-stream  wind  components  had  zero  means,  and 
wrote  the  data  on  magnetic  tape  for  further  analysis.  The  mean,  rms, 
maximum  and  minimum  values  and  the  trends  were  printed,  as  well  as  the 
momentum  and  heat  fluxes  and  the  drag  coefficient  for  a wind  extrapolated 
to  a height  of  10  m above  the  surface.  The  extrapolation,  which  assumes 
wind  velocities  vary  logarithmically  with  height,  introduced  small  corrections 
and  negligible  errors  in  the  10  m wind  velocity  values,  U . 
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A subroutine  (GILL)  was  added  to  process  the  data  from  the  UBC  Gill 
three-component  anemometer,  and  a small  number  of  associated  changes  were 
made  in  the  main  program. 

Our  spectral  analysis  programs  (Dobson  et  al. , 197*+)  read  data 
from  program  A TO  D,  which  removes  means  and  trends.  Detrending  removes 
some  of  the  variance  associated  with  periods  longer  than  the  run  length. 

This  slightly  reduces  the  downwind  turbulence  levels  which  we  report,  but 
the  part  which  is  removed  was  not  adequately  sampled  in  the  first  place. 

The  effect  on  vertical  winds  and  hence  the  turbulent  fluxes  is  negligible. 

The  highest  frequency  analysed,  5 Hz,  is  determined  by  the  sampling  rate 
of  10  per  second.  The  lowest  frequency,  0.005  Hz  (sometimes  0.01  Hz), 
is  determined  by  the  length  of  blocks  of  data  used  in  the  spectral  analysis. 

Dubious  data  have  been  deleted  from  the  results  presented  here 
(Table  l).  The  sonic  anemometer  tended  to  mistrigger  in  damp  weather 
causing  'spikes'  in  the  output  which  led  to  rejection  of  a number  of  data 
runs.  The  Mark  7 thrust  anemometer  suffered  a connector  failure,  thus 
missing  the  first  nine  runs.  Runs  16  and  35  were  not  digitized  because 
they  were  recorded  when  the  wind  was  blowing  down  off  the  dunes . Other 
reasons  for  deletions  included  excessive  vibration  of  the  thrust  anemometer 
or  signals  running  off-scale. 

4.  RESULTS 

4.1  Velocity 

The  sonic  anemometer  will  be  the  basis  for  comparison  of  results 
because  it  is  more  widely  used,  has  a linear  response  which  permits  straight- 
forward data  processing,  and  has  a flat  response  over  our  frequency  range 

0 to  5 Hz. 
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Figure  3 shows  the  mean  wind  speeds  from  the  Gill  anemometer  (U  ) 

uILL 

UTHRUST  T>1  and  UTHRUST  6'U  plotted  versus  US0NIC.  A linear  re8ression  for 
the  average  wind  speed  at  10  m,  U , gives 


(U10^GILL  = 1,02  (U10^S0NIC  + 0,33  m/s' 


The  sonic  anemometer  gives  lower  values  for  the  mean  wind  speed 
than  the  Gill  anemometer  (Fig.  3).  Calibrations  of  the  sonic  anemometer 
electronics  by  introducing  measured  delays  in  its  sound  pulses  (Table  2) 
showed  it  to  be  within  the  manufacturer's  + 3%  specification  of  the  factory 
values  (3,  10,  and  30  m/s  full-scale  for  horizontal  components  and  1,  3, 

10  m/s  for  vertical  components)  and  so  didn't  change  by  more  than  2%  between 
calibrations . 

The  higher  wind  speeds  indicated  by  the  UBC  Gill  anemometer  may  be 
attributed  to  small  errors  in  calibrating  the  anemometers  and  in  setting 
the  voltage  corresponding  to  zero  wind  speed. 

The  Gill  anemometer  is  calibrated  by  rotating  the  propeller  shaft  at 
known  rates  and  measuring  the  output  voltage.  The  conversion  factor 
between  rotation  rate  and  wind  velocity  is  supplied  by  the  factory  and, 
although  not  done  for  this  particular  anemometer,  direct  wind  tunnel 
calibrations  of  Gill  anemometers  at  UBC  and  elsewhere  have  always  shown 
excellent  agreement  with  factory  calibrations  (l-2$).  Post-experiment 
calibrations  confirmed  the  earlier  ones  within  experimental  error  (about 
one  part  in  250).  Uncertainty  in  the  dc  level  at  zero  wind  and  in  the 
offset  are  equivalent  to  0.005  and  0.15  m/s  respectively.  The  first 
could  be  systematic,  while  0.1  m/s  of  the  second  is  expected  to  be 
random.  Thus,  the  Gill  anemometer  ought  to  indicate  velocities  that  are 


at  most  2%  high. 
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TABLE  2 

SONIC  ANEMOMETER  NO.  2 CALIBRATION 

(FULL  SCALE  RANGE  AT  20C) 
PERFORMED  BY  TIMING  PULSES 


Rarif'e 

(Be fore 

(09/09/75)* 

After 

(18/11/75) 

A 

B 

W 

A 

B 

W 

m/s 

m/s 

m/ s 

m/ s 

m/s 

m/s 

1 

3.03 

3.00 

1.00 

2.99 

3.01 

0.98 

O 

(— 

9.80 

9.91 

3. Oit 

9-jk 

9.87 

3.02 

3 

29.07 

30.28 

9.89 

28.92 

29.50 

9.80 

* These  values  used  in  data  analysis. 
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The  response  of  the  Gill  anemometer  propellers  falls  off  at  high 
wave  numbers,  and  is  3 dB  down  at  a fluctuation  wavelength  of  about  6 m. 

The  turbulence  scales  (wavelengths)  that  contribute  to  the  variances  and 
covariances  are  proportional  to  height  and  at  a fixed  height  are  very 
nearly  independent  oi  wind  speed.  The  reduction  of  measured  variances  and 
covariances  due  to  anemometer  response  limitations  is  thus,  fortunately, 
also  independent  of  wind  speed.  At  the  height  of  observation,  about  10  m, 
we  expect  the  indicated  covariance  u^u^  to  be  about  8s%  of  its  true  value. 
Examination  of  several  cospectra  showed  that  for  the  Gill  about  15% 

of  the  stress  was  lost  due  to  small  scale  response  limitations.  Also,  the 
stress  from  the  Gill  and  Sonic  agree  on  average  if  this  condition  is  made 
and  the  calibration  difference  between  the  two  instruments  (discussed  later) 
is  taken  into  account.  Since  the  BIO  data  processing  system  does  not  lend 
itself  to  making  frequency-dependent  corrections,  we  report  indicated 
values  for  velocity  fluctuation  components  keeping  in  mind  that  the  vertical 
wind  variance  will  be  reduced.  Corrected  values  of  the  covariance  u^u^  are 
used  in  Table  1.  While  errors  from  either  anemometer  alone  cannot 
account,  for  the  difference  between  the  sonic  and  the  Gill  anemometers, 
the  combined  possible  errors  of  the  two  anemometers  can. 

Thrust  anemometer  calibrations  (Table  3)  made  before  the  experiment 
using  the  200  ft  field  cable,  and  used  in  the  analysis,  indicated  larger 
full-scale  ranges  (lower  sensitivity)  than  calibrations  made  after  the 
experiment  without  the  long  cable.  The  post-experiment  calibration  of 
Thrust  G.h  agrees  within  1%  in  velocity  with  each  of  a series  of  seven 
calibrations  from  March  to  August  197^.  The  nnges  of  wind  speeds  observed 
with  the  thrust  anemometers  were  not  sufficient  for  regression  analysis  but 
the  observed  wind  speeds  agreed  better  with  the  Gill  anemometer  than  with 


the  sonic  anemometer: 


THRUST  AUDIOMETER  CALIBRATION 
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r- 


UTHRUST  6.l/UGILL  ' 
UTHRUST  T.1/UGILL 
UTHRUST  6.1+^US0NIC 
UTHRUST  T.l^SONIC 


0.990  ±0.027  standard  deviation 
0.984  ±0.016 
= 1.042  ±0.009 
= 1.028  ±0.015 


(12  runs) 
( 6 runs ) 
( 5 runs ) 
( 5 runs ) 


These  results  suggest  that  the  Gill  anemometer  was  reading  slightly  high 
(l  to  2%)  while  the  sonic  anemometer  was  slightly  low  (3  to  h%) . Although 
the  thrust  anemometers  are  not  necessarily  more  accurate  than  the  others  they 
were  calibrated  in  a wind  tunnel  while  the  others  were  calibrated  less 
directly,  by  timing  of  pulses  (sonic)  and  by  rotation  rate  (Gill  anemometer). 
In  preparation  for  the  work  on  the  stable  platform,  calibration  of  the  Gill, 
thrust  and  Bendix  Aerovane  are  being  done  in  the  wind  tunnel  at  BIO  to  try 
to  reduce  calibration  differences. 

Comparing  the  standard  deviation  of  downstream  velocity  (Fig.  4)  we 
see  that  the  sonic  anemometer  generally  indicates  slightly  lower  levels 
than  the  Gill  anemometer,  by  about  the  same  amount  as  the  difference  in 
indicated  mean  wind. 

The  lateral  response  of  the  Gill  anemometer  is  determined  by  the 
characteristics  of  the  vane,  which  apparently  tends  to  overshoot  and 
increase  the  indicated  turbulence  levels  by  a further  5$  (Fig.  5).  Thrust 
anemometer  7.1  suffers  badly  from  lateral  vibration,  while  even  thrust  6.4 
which  is  counterbalanced  against  horizontal  vibration  gives  values 
some  15  to  20%  higher  than  those  from  the  sonic  anemometer.  The  guy  wires 
supporting  the  instrument  mast  were  observed  to  sing  in  high  winds  and  this 


appears  to  have  had  a detrimental  effect  on  the  thrust  anemometers. 
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Vertical  turbulence  indicated  by  the  Gill  anemometer  is  only 
about  75$  of  that  indicated  by  the  sonic  anemometer  (Fig.  6),  as  a result 
of  attenuation  of  the  smaller  scales.  Mounting  the  Gill  anemometer  on  a 
higher  tower,  as  was  originally  planned,  would  reduce  this  problem.  As 
noted  before  the  reduction  is  essentially  independent  of  wind  speed  so  a 
correction  is  easily  made.  Both  thrust  anemometers  tend  to  give  higher 
a ^ values  than  the  sonic  anemometer  by  an  amount  comparable  to  the 
difference  in  indicated  mean  wind. 

To  test  the  ability  of  the  anemometers  to  resolve  whether  the  wind 
is  relatively  more  or  less  gusty,  we  may  compare  the  indicated  turbulence 
levels  divided  by  the  mean  wind  speed.  A linear  regression  from  data  in 
Table  1 gives  very  good  agreement  for  the  downwind  compon  nt . 


'WoiIA  = °-003  * °'97  'WsOKIC  S 0M- 


correlation  coefficient  0.96 


The  distribution  of  points  and  the  above  regression  line  are  not  the  same 
as  in  Figure  b,  since  runs  at  higher  wind  speeds  must  have  higher  o values 
but  need  not  have  higher  o^/U  values.  The  thrust  anemometers  indicated 
higher  levels  with  the  mean  ratio  ( a i /U)THRUST//  1 /^SONIC  being 

1.10  to.oH  and  1.13  ±0.10  for  Thrust  G.b  and  Thrust  7-1  respectively. 

In  the  lateral  turbulence  levels  a similar  pattern  emerges,  but  the 
agreement  is  not  quite  so  close,  a regression  line  being 

<°2/U)CILL  ' °‘019  * °'73  <“2/UWc 


correlation  coefficient  0.91 


O],  ( m /s  ) 


j 


SONIC  0£(m/s) 


Figure  5-  Comparison  of  cross-wind  turbulence  levels,  with  regression 

line  °2  GILL  = -9*013  m/s  + X-lU3  °2  SONIC’ 
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Thrust  7.1  suffered  particularly  badly  from  lateral  vibration  with  the 
ratio  (°o/u)thrust/(02/U)sONIC  = 1,T6  ±0-35  while  for  Thrust  6.1  the  ratio 
was  a more  acceptable  1.21  ±0.13.  The  lateral  turbulence  is  not  required 
in  our  eddy  flux  calculations. 

A linear  regression  comparison  of  vertical  turbulence  levels  from 
Table  1 gives 

(03/UW  ' 0'011  + °'1*8  (°3/U)S0»IC  10-002 
correlation  coefficient  0.34 

The  first  term  (intercept)  accounts  for  about  one-quarter  and  the  second 
term  for  about  one-half  of  the  rms  vertical  velocity  indicated  by  the  sonic 
anemometer,  the  remaining  one-quarter  being  lost  because  of  small-scale 
response  limitation  of  the  Gill  anemometer.  The  correlation  between  Gill 
and  sonic  anemometers  is  not  nearly  so  good  for  the  vertical  as  for  the 
downwind  component,  but  this  can  be  explained  by  the  fact  that  o^/U  varies 
less  than  the  horizontal  turbulence  levels  so  that  relatively  small  sensor 

G 

variations  become  more  apparent.  The  thrust  anemometers  agreed  well  with 
the  sonic  anemometer  with  ratios  (CT3/U)TrousT/(a3/U)s0NIC  of  1*02  ±0-°9  and 
0.99  ±0.07  for  Thrust  6.4  and  Thrust  7*1  respectively. 

4.2  The  brag  Coefficient 

The  drag  coefficient,  = -u^u  /(U.^)  , where  u^  and  u^  are  the 
downwind  and  vertical  wind  velocity  fluctuations  respectively,  has  been 
calculated  for  the  four  anemometers  and  is  plotted  in  Figure  7.  In  this 
plot,  values  used  in  determining  each  were  obtained  from  the  same 
instrument,  so  the  effects  of  calibration  differences  between  the  instru- 
ments should  be  reduced.  As  noted  earlier  the  values  of  u^u^  reported 
for  the  Gill  have  been  corrected  for  its  small  scale  response  limitations 
by  using  originally  calculated  value  divided  by  O.85. 

. _ j 


4 
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Although  there  is  scatter,  the  results  agree  well  on  the  whole. 

Since  a possible  error  of  1°  in  the  orientation  of  the  coordinate  system 

produces  an  error  of  10$  in  the  value  of  u^u^,  for  any  one  run  values  from 

two  systems  could  easily  differ  by  20$.  Detrending  before  calculating  u^ 

can  also  affect  the  u^  and  u^u^  values  from  the  Gill  anemometer.  It  appears 

that  errors  introduced  by  detrending  when  the  trend  and  variation  are  large 

may  cause  low-frequency  contributions  to  u_,u_  nTTT  to  be  too  small. 

J.  5 ulLL 

Least  squares  linear  fits  to  the  data  for  the  drag  coefficients 
from  the  sonic  and  Gill  anemometers  give 


(103C 

(103C 


10  ^ SONIC 
10 ^ GILL 


= 0.59  + 0.063  (U 
= 0.1+1  + 0.083  (U 


10  ^ SONIC 
10 ^ GILL 


Because  of  the  few  points  and  small  range  of  wind  speed,  regression  line 
fits  for  the  thrust  anemometers  were  not  warranted;  however,  as  can  be  seen 
in  Figure  7»  the  results  are  in  reasonable  agreement  with  those  from  the 
sonic  and  Gill  anemometers.  The  overall  result  is  consistent,  considering 
scatter,  with  the  relation  103C1Q  = 0.63  + 0.066  U reported  by  Smith  and 
Banke  (1975)* 

The  values  for  C averaged  over  all  wind  speeds,  shown  in  Table  1, 

are  not  inconsistent  with  one  another  and  with  previous  data  (e.g.  Smith 

and  Banke,  1975)-  Much  of  the  difference  in  values  is  associated  with 

averaging  over  different  runs.  For  20  runs  for  which  the  Gill  and  sonic 

3 

anemometers  both  give  values,  10  = 1.19  ±0.24  and  1.18  ±f\022  (mean  ± 

standard  deviation)  for  Gill  and  sonic  anemometers  respectively.  For  17 
runs  with  Gill  and  thrust  anemometers,  the  values  are  respectively  1.37 
±0.25  and  1.49  ±0.26.  For  11  runs  with  sonic  and  either  of  the  thrust 
anemometers  the  values  are  1.23  ±0.24  and  1.4l  ±0.27  respectively.  The 
thrust  anemometer  values  are  somewhat  high,  perhaps  in  part  due  to  the 


a SONIC 
O THRUST  6.4 
• THRUST  7.1 


U,0(  m/s ) 


Figure  7-  Drag  coefficients  plotted  against  wind  speed  with  regression 

lines  for  Gill  and  Sonic  and  the  regression  line  from 
Smith  and  Banke  (1975). 
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vibration  noise.  Very  high  thrust  anemometer  drag  coefficients  for  runs  18 
and  33  were  not  used  in  calculating  means  and  standard  deviation. 

The  fact  that  not  including  a few  runs  can  change  C by  10  to  20$ 
even  when  using  the  same  site  and  instrument  is  worth  noting.  Much  of  the 
scatter  in  existing  results  and  differences  among  them  may  be  associated 
with  having  small  samples  of  a variable  process . 

The  values  for  o_/u»  and  the  overall  correlation  coefficients 

3 * 

..... 


-r^  = -u^u^/Co^o^) 2 given  in  Table  1 are  within  the  usual  range  of  values 

for  turbulence  in  the  lower  boundary  layer.  Smith  and  Banke  (1975)>  for 

example,  found  o_/u„  = 1.2  ±0.1  and  r. _ = -0.3U  ±0.07.  The  Gill  anemometer 
3 * 13 

values  are  slightly  different  from  the  rest.  Comparison  with  o 

j oUIM_LU 

and  examination  of  <t>^  spectra  show  that  the  Gill  underestimates  by 
25  to  30$  because  of  its  small  scale  response  limitations.  If  we  take 
the  correct  value  to  be  1.35  cr_  r,TTT  then  the  average  a_/u.  becomes  1.31 

3 uiliii  3 

and  the  average  -r  becomes  0.29,  in  much  better  agreement  with  values 


from  the  sonic  and  thrust  anemometers.  Note  that  in  Table  1 u^u^ 
has  been  corrected  for  response  limitations  but  has  not. 

4 . 3 iipectra 

Spectra  and  cospectra  of  the  velocity,  temperature,  and  humidity 
fluctuations  have  been  analyzed  over  a frequency  band  from  0.005  Hz 
(sometimes  0.01  Hz)  to  5 Hz.  Logarithmically  spaced  bandwidths  have  been 
chosen.  Examination  of  the  spectra  can  reveal  deficiencies  of  the 


sensors  such  as  noisy  operation  or  limited  frequency  response,  since  the 
general  shapes  of  the  spectra  are  well  known  from  previous  studies. 

Velocity  spectra  for  Run  32  are  typical  of  data  obtained  with 


all  three  anemometers  functioning  well.  The  log-log  plots  (Fig.  8a,c,e) 
are  useful  for  showing  power-law  spectral  shapes  as  straight  lines. 
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while  plots  of  (frequency  x spectrum)  against  logarithmic  frequency 

(Fig.  8b,d,f)  have  areas  proportional  to  kinetic  energy  atul  so  give  a 

better  Indication  of  the  relative  importance  of  spectral  features.  The 

downwind  velocity  components  (Fig.  8a, b)  show  reasonable  agreement,  with 

the  UBC  Gill  anemometer  spectrum  lying  between  the  other  two  except 

above  0.5  Hz  where  its  frequency  response  falls  off.  Different  frequency 

bands  were  chosen  for  the  thrust  anemometer,  giving  slightly  different 

relative  sizes  of  small  peaks  within  the  spectrum.  The  expected  -5/3 

slope  is  seen  in  the  log-log  plot  from  0.2  Hz  to  2 Hz  in  the  thrust  and 

sonic  anemometer  spectra,  while  'aliasing'  of  higher- frequency  signals  increases 

the  computed  spectra  at  frequencies  around  3 to  5 Hz.  The  Gill  anemometer 

spectrum  has  a slope  of  -5/3  over  a relatively  small  range  from  0.2  to 

0.5  Hz. 

The  lateral  or  crosswind  spectra  (Fig.  8c, d)  show  better  agreement 
in  small  peaks  and  valleys  since  all  three  were  analyzed  with  the  same  array 
of  bandwidths.  The  UBC  Gill  anemometer  spectrum  is  higher  than  the 

others  in  the  range  0.1  to  1 Hz,  possibly  due  to  overshooting  of  its  direction 

vane,  while  its  response  falls  off  rapidly  above  1 Hz.  The  thrust  anemometer 
gives  higher  values  than  the  sonic  at  all  frequencies.  All  the  spectra  have 
a very  small  peak  near  1 Hz  which  may  be  due  to  vibration  of  the  tower. 

In  the  vertical  direction  (Fig.  8e,f)  the  response  of  the  UBC  Gill 
anemometer  is  3 dB  down  at  a frequency  of  about  1 Hz,  accounting  for  the 

lower  values  of  mentioned  earlier.  The  sonic  anemometer  spectrum  shows 

a -5/3  power  law  for  a small  range  of  frequencies  around  1 Hz  but  this 
behavior  is  not  nearly  as  well  defined  as  for  the  downwind  component.  The 
range  of  frequencies  analyzed  barely  reaches  the  area  where  local  isotropy 
and  the  associated  f-5/3  power  law  behavior  is  possible.  The  observed 
spectra  generally  follow  the  expected  shape.  At  1 Hz  the  ratios  4>^^(  f)/4>1^(f) 


d stress  cospectj 
ic  anemometer  ( si 
e),  and  thrust  anemometer  (dotted  line). 
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for  the  sonic  and  thrust  anemometers  are  1.36  and  l.lU  respectively  for  this 


r 

I 

I 


i 


run.  Local  isotropy  theory  predicts  that  this  ratio  should  approach  U/3  at 
the  higher  frequencies. 

The  wind  stress  cospectra  ftfi^lf)  (Fig.  9)  show  general  agreement 
in  siiape,  with  the  thrust  anemometer  indicating  higher  stress  than  the 
other  two.  Small  values  at  the  low  and  high  frequency  end  of  the  cospectrum 
indicate  that  the  frequency  band  analyzed  contains  most  of  the  range  of 
interest  for  the  wind  stress  measurement. 

Coherence  and  phase  between  measurements  of  the  velocity  components 
show  good  agreement  at  the  lower  frequencies , but  at  higher  frequencies 
there  is  a rapid  drop  in  coherence  accompanied  by  a more  gradual  shift  in 
phase.  In  Run  32  (Fig.  10)  the  thrust-sonic  anemometer  coherence  persists 
to  higher  frequencies  partly  because  the  separation  between  these  two  ane- 
mometers (l  m versus  1.7  m)  is  less  than  that  of  the  UBC  Gill  and  sonic 
anemometers.  For  this  run  the  sonic  and  thrust  anemometers  and  the  UBC  Gill 
propellers  were  separated  in  a nearly  crosswind  direction.  Nearly  all  of  the 
observed  phase  shift  for  u^  and  u^  is  due  to  the  response  characteristics  of 
the  Gill  anemometer.  The  vane  from  which  u0  ~TTT  is  calculated  is  almost 
1 m downstream  relative  to  the  sonic  and  so  of  the  observed  u0  phase  dif- 
ference (90°  at  1 Hz),  nearly  half  (U0°  at  1 Hz)  may  be  due  to  this  displace- 
ment. The  thrust  and  sonic  anemometers  remain  in  phase  (±10°)  up  to  frequencies 
at  which  the  coherence  drops  to  small  values.  The  half-coherence  frequency 
between  thrust  and  sonic  anemometers  is  lower  for  the  downwind  component 
(0.8  Hz)  than  it  is  for  the  lateral  and  vertical  components  (1.5  Hz). 

Because  the  data  were  sampled  10  times  per  second,  but  low-pass  filtered 
only  by  3 db  at  11  Hz  (18  dB  per  octave  rolloff),  energy  at  frequencies  above 
5 Hz  is  folded  back  of  "aliased"  into  the  5 Hz  band  covered  by  the  spectral 
analysis.  The  thrust  anemometers  were  particularly  noisy  at  high  frequencies, 
and  aliasing  may  explain  their  generally  higher  spectral  levels. 


Coherence  spectra  for  Run 
32  between  UBC  Gill  ane- 
mometer and  sonic  anemom- 
eter (dashed  line)  and 
between  thrust  anemometer 
and  sonic  anemometer  (dotted 
line).  Phase,  relative  to  sonic 
anemometer-,  of  UBC  Gill  anemometer 
(X)  and  of  thrust  anemometer  (0). 
(a)  downwind,  (b)  cross-wind, 
and  (c)  vertical. 
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It . It  Temperature  Fluctuations  and  Sensible  Heat  Flux 

The  mercury  thermometer  readings  for  air  temperature,  T,  (Table  4 ) 
were  taken  near  the  top  of  the  tower  either  before  or  after  recordings . 

Thus  in  addition  to  reading  error  of  ±0.1  to  0.2C  there  is  uncertainty 
because  they  are  not  averages  over  the  run  as  are  the  microbead  values. 

Since  o is  of  the  order  of  0 . 1 | AT | , the  difference  between  instantaneous 
and  mean  temperatures  could  easily  be  ±0.2  AT.  This  estimate  is  based  on 
the  assumption  of  Gaussian  statistics  and  no  large  trends.  Some  of  the 
differences  between  the  mercury  thermometer  and  the  microbead  values  could 
be  due  to  these  errors. 

Allowing  for  ±0.5%  of  full-scale  nonlinearity  or  drift,  the  BIO 
thermistor  should  have  an  accuracy  of  ±0.4C.  The  mean  temperature  from 
the  BIO  thermistor  agreed  within  1C  with  the  mercury  thermometer  readings 
(Table  4)  except  for  runs  28,  29  and  30  where  the  BIO  thermistor  behaved 
erratically,  as  will  be  discussed  later. 

When  it  is  working  properly  the  UBC  thermistor  should  give  an 
average  air  temperature  to  ±0.2C.  The  UBC  thermistor  and  the  mercury 
thermometer  agreed  within  0.3C  for  the  first  few  runs  (2,  6,  7»  8 and  9 
in  Table  4),  while  afterward  the  differences  were  generally  larger  than 
1C.  If  the  mean  is  not  correct  the  fluctuations  may  also  have  some  errors. 
We  suspect  shunting  of  the  bead  by  salt  films  and  in  future  shall  insulate 
the  leads  and  have  a glass  rod  thermistor  near  the  microbead  to  check  it. 

A water  sample  was  taken  from  the  surf  with  a bucket  and  T obtained 
from  it  with  a reading  error  of  ±0.2C.  Thus  the  errors  in  T are  perhaps 
±o.3C  from  reading  plus  20%  from  sampling  and  are  quite  large. 

Comparing  the  o^/AT  values  from  the  two  systems,  we  see  that  the  UBC 
values  are  smaller.  Negative  values  are  for  cases  with  negative  At.  When 
|AT|<0.5C  (runs  9»  11,  12,  23,  28  and  31)  the  standard  deviation  of  temper- 
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TABLE 


1+ 


TEMPERATURE  DATA 


Run 

Thermometer 

Sonic 

Gill 

BIO 

Thermistor 

UBC 

Thermistor 

T 

Ts 

uat 

uat 

T 

°t 

ugt 
C m/s 

T 

°t 

u^t 

C 

C 

C m/s 

C m/s 

C 

AT 

C 

AT 

C m/s 

(or  0) 

(or  0 ) 

t 

1 

2 

18.0 

16.6 

- 

-3.9 

17.2 

-0.il* 

-0.009 

18.0 

-0.12 

0.000 

6 

13.5 

lU.O 

3.8 

- 

12.9 

0.21+ 

0.013 

13.3 

0.11+ 

- 

T 

13.2 

lU.O 

5.8 

6.1+ 

13.0 

0.51 

0.012 

13.1* 

0.10 

0.005 

0 

13.3 

lU.O 

5.9 

6.1* 

13.2 

0.20 

0.019 

13.6 

0.11 

0.008 

9 

13.7 

13.8 

0.8 

0.8 

13.6 

(0.16) 

0.013 

ll+.O 

(0.11) 

0.006 

10 

12.5 

11.5 

-9.8 

-10.3 

11.8 

-0.08 

0.009 

13.8 

-0.08 

0.006 

11 

12.3 

12.2 

-1.1 

-1.2 

11.9 

(0.08) 

0.012 

ll+.l* 

(0.08) 

0.006 

12 

12.6 

12.3 

-3.1 

-3.3 

12.1 

(0.09) 

0.007 

15.1 

(0.07) 

O.OOL 

13 

12.8 

12.3 

-5.5 

-5.8 

12.6 

-0.12 

0.007 

15.2 

-0.06 

0.003 

14 

13.0 

12.2 

- 

-10.1 

12.5 

-0.08 

0.007 

15.1* 

-0.06 

0.005 

15 

12.6 

12.1 

- 

-6.1* 

12.6 

-0.12 

0.007 

15.8 

-0.07 

0.003 

18 

9.6 

10.7 

- 

13.8 

9.5 

0.23 

0.01*6 

8.3 

0.15 

0.029 

19 

10.0 

11.1 

- 

13.6 

9.8 

0.32 

0.01+3 

8.7 

O.ll* 

0.021+ 

20 

10.0 

11.2 

- 

lit. 8 

9.1* 

0.26 

0.031+ 

9.0 

0.11 

0.020 

21 

10.8 

11.3 

6.0 

6.1* 

9.8 

0. 71+ 

0.026 

9.1* 

0.22 

0.018 

23 

11.2 

11.0 

- 

-2.9 

10.8 

(0.30) 

0.009 

9.8 

(0.09) 

0.008 

2h 

11.2 

11.0 

- 

-3.0 

- 

- 

- 

9.3 

(0.06) 

0.008 

20 

12.0 

11.6 

- 

-1+.8 

10.1* 

(0.10) 

- 

9.8 

(0.05) 

0.001 

29 

12.3 

12.0 

-2.9 

-3.0 

11.1 

- 

-0.006 

10.2 

(0.01+) 

-0.001 

30 

12.6 

12.0 

-5.6 

-5.7 

11.3 

-0.10 

-0.002 

10.3 

-0.05 

-0.001 

31 

12.6 

12.0 

-5.3 

-5.1* 

11.7 

-0.06 

-0.003 

10.1* 

-0.01+ 

-0.001 

33 

12.7 

12.2 

-5.5 

-5.8 

13.0 

-0.10 

-0.005 

10.1* 

-0.06 

-0.002 

Note : 

ot(C) 

is  given  in 

brackets 

in  place 

of  o^/AT 

in  cases  where 

1 AT | <0 . 

■ 5C 
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ature  is  listed  in  Table  4 with  brackets  to  allow  comparison  of  values 
from  the  two  systems.  The  0 / aT  values  are  not  nearly  so  uniform  as  those 

Xt 

for  wind  components.  Much  of  this  scatter  could  be  associated  with  the 
errors  in  AT.  Because  of  the  small  aT's  the  signals  are  small  and  there 
may  be  some  noise  included.  Finally,  there  are  some  cases  where  one  or  the 
other  system  is  not  working  properly.  For  example,  the  BIO  thermistor  data 
for  Run  21  are  dubious  since  large  indicated  fluctuations  of  3-  to  5-minute 
periods  were  not  corroborated  by  the  UBC  thermistor.  These  fluctuations  showed 
up  as  a large  value  ofo^/AT  (0.74),  high  spectral  values  at  low  frequencies, 
and  reduced  coherence  at  low  frequencies  between  the  BIO  and  UBC  uhermistors. 

This  phenomenon  persisted  through  Run  24,  followed  by  excessive  drift  of  the 
BIO  thermistor  during  Runs  25  and  26.  Rain  fell  during  Runs  23  to  26,  and 
a possible  explanation  of  this  peculiar  behavior  is  salt  or  moisture  contam- 
ination of  the  BIO  thermistor  bead. 

Even  though  the  microbeads  are  similar,  an  examination  of  the  temper- 
ature spectra  (Fig.  11a)  shows  higher  values  for  the  BIO  thermistor  at  all 
frequencies.  At  high  frequencies  the  values  of  the  UBC  thermistor  fall  off 
more  rapidly  and  its  phase  lags  the  BIO  thermistor  (Fig.  lib).  This  is 

; ti 

probably  due  to  the  shield  which  kept  the  UBC  sensors  out  of  the  direct 
airflow.  The  sensor  location  has  been  changed  to  improve  this  response 
in  future  measurements.  The  shield  may  also  contribute  to  the  general 
lower  and  o / AT  values  of  the  UBC  thermistors.  However,  it  may  be 
that  the  BIO  thermistor  values,  particularly  the  peak  at  .03  Hz,  are  too 
high.  A possible  cause  is  occasional  small  step  changes  in  the  servo- 
potentiometer  which  balances  out  the  mean  temperature.  As  previously 
mentioned,  salt  or  moisture  contamination  of  one  or  both  of  the  thermistor 
beads  is  always  possible  and  could  also  have  contributed  to  these  effects. 

Error  bars  in  Fig.  11a  give  90?  confidence  intervals  for  the  spectral 


estimates  and  are  generally  similar  to  the  error  bars  (not  shown)  for  spectra 
of  other  variables  in  this  paper. 
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Even  if  the  values  are  affected  by  noise  it  may  still  be  pos- 
sible to  measure  the  covariance  u^t  (e.g.  Run  23,  Table  i+)  since  the  noise 
in  the  two  channels  probably  has  a low  correlation.  Some  of  the  recording 
noise,  that  associated  with  tape  speed  variation,  will  be  correlated;  compen- 
sation is  done  to  reduce  this  noise  but  it  cannot  be  eliminated  entirely. 

Figure  12a  shows  u^t  from  the  two  systems  plotted  against  one 

another.  These  values,  taken  from  Table  4,  have  not  been  corrected  for 

sensor  response  limitations.  Most  of  the  values  are  fairly  small  but  there 

does  seem  to  be  a consistent  relationship  particularly  for  the  larger  values. 

It  appears  that  (u_t)  Trnri  * 0.6  (u  t)  . Much  of  this  difference  is 
j UBL  j nlu 

probably  due  to  the  response  limits  of  the  UBC  sensors  although  some  of  it 
may  also  be  due  to  calibration  errors  and  noise. 

Figure  12b  shows  u^t  from  each  system  plotted  against  UAT  with  U 
from  the  same  system  and  AT  from  the  mercury  thermometer  observations  as 
noted  before.  There  is  a great  deal  of  scatter  which  is  caused  by  such 
things  as  the  large  possible  errors  in  AT  values  and  the  small  values  of 

u.^t  leading  to  possible  noise  problems.  The  slopes  of  the  regression  lines 

-3  -3 

are  1.3  x 10  and  1.0  x 10  for  the  BIO  and  UBC  data  sets,  respectively. 
Most  of  this  difference  is  probably  due  to  response  limitations  noted 

before;  a correction  based  on  the  slope  of  Figure  12a  would  give  a slope 

-3 

of  1.7  x 10  for  the  UBC  data.  Both  sets  of  data  show  an  upward  flux 
u.,t  = 0.0095  Cm/s  (BIO)  and  0.0076  Cm/s  (UBC)  at  UAT  = 0 which  may  be 
caused  by  the  scatter,  AT  errors,  noise,  or  may  be  real. 

The  slopes  of  the  regressions  are  consistent  with  recent  observa- 
tions of  C^,  = [u^t/ (UAT) ] reported  by  Friehe  and  Schmitt  (1976)  although, 
in  view  of  the  small  values  of  u^t  and  the  scatter,  our  values  cannot  be 


r 


considered  very  reliable.  However,  the  temperature  fluctuations  are  small 


I03  u3t  ( BIO  ) °C  m/s 


Figure  12.  (a)  Comparison  of  u t from  the  BIO  and  UBC  systems.  Th 

line  has  a slope  of  0.6. 

(b)  Heat  flux  u^t  as  a Function  of  UAT,  with  regression 
lines  for  BIO  thermistor  data  (solid)  and  UBC  data 
(dashed) . 
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and  the  sensible  heat  flux  is  much  smaller  than  the  latent  heat  flux 
(when  it  was  measured,  and  hence  much  of  the  time  since  the  overall  conditions 
were  similar).  Under  such  conditions  the  temperature  signals  may  be  con- 
taminated by  the  response  of  saline  films  on  the  thermistors  to  moisture 
fluctuations  (Friehe  and  Schmitt,  1976).  This  would  give  u^t  values  that 
are  too  large,  and  apparent  cospectra  that  have  relatively  larger  values  at 
higher  frequencies  than  do  the  u^u^  cospectra.  Such  effects  could  explain 
the  positive  intercepts  in  Fig.  12b  and  the  scattered  values. 

4 . 5 Humidity 

During  a four-hour  period,  satisfactory  recordings  of  Lyman-alpha, 
Brady  and  aluminum  oxide  humidity  data  were  made  (Runs  6 to  9,  Tables  1 and 
5).  Unfortunately,  there  was  an  offset  change  for  the  Gill  velocity  values 
near  the  start  of  Run  6 and  the  quantities  dependent  on  velocity  cannot  be 
used  for  this  run  as  analyzed.  Reasonable  correlations  of  mean  humidity 
between  the  three  sensors  were  obtained,  but  the  sensors  had  some  calibra- 
tion drifts  and  the  UBC  sensors  had  limited  frequency  response. 

The  sensors  were  separated  90  cm  horizontally  and  85  cm  vertically. 
Sonic  u^  was  used  with  the  Lyman-alpha  data  while  Gill  u^  was  used  with  the 
Brady  and  aluminum  oxide  data  to  derive  the  moisture  flux,  qu^.  The  data 
in  Table  5 are  not  corrected  for  response  limitations. 


The  spectrum  from  the  unshielded  Lyman-alpha  (ultraviolet)  absorp- 
tion humidiometer  for  Run  6 (Fig.  13a)  resembles  that  of  downwind  velocity 
and  shows  an  approximate  -5/3  power  law  starting  at  0.3  Hz.  The  two  UBC 
humidity  sensors  showed  much  lower  humidity  fluctuation  levels  at  all  fre- 
quencies and  also  additional  attenuation  with  frequencies  from  0.02  Hz  up. 
Of  the  two,  the  Brady  type  sensor  had  the  higher  spectral  levels  at  0.01 
to  0.02  Hz,  but  its  response  fell  off  more  rapidly  at  frequencies  from 


0.02  to  0.1  Hz. 


f (Hz) 

Figure  13.  (a)  Humidity  spectra  for  Run  6;  Lyman-alpha  (solid); 

Brady  (dashed),  and  AlO^  sensor  (dotted). 

(b)  Coherence  between  Brady  and  Lyman-alpha  (dashed  line) 
and  between  A1C>2  and  Lyman-alpha  (dotted  line).  Phas 
relative  to  Lyman-alpha,  of  Brady  (X)  and  A1C>2  (0) 
humidity  signals. 
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The  coherence  between  the  Lyman-alpha  and  the  other  sensors  is 
good  up  to  0.17  Hz  but  a phase  lag  of  90°  develops  by  0.13  Hz  in  the  Brady 
and  0.3  Hz  in  the  aluminum  oxide  sensor  (Fig.  13b).  Very  little  of  this 
phase  shift  can  be  accounted  for  by  sensor  separation.  Thorpe  (1973) 
found  that  there  is  a delay  other  than  a simple  R-C  filter  type  smoothing 
of  the  fluctuations  measured  by  a Brady  relative  humidity  sensor  and 
that  this  could  not  be  accounted  for  by  temperature  effects  in  the  sensor. 

Neither  of  the  UBC  humidity  sensors  in  their  present  form  have 
adequate  frequency  response  for  eddy  flux  measurements  of  the  type  reported 
here.  The  aluminum  oxide  sensor  had  large  calibration  shifts  because  of 
salt  contamination  and  will  no  longer  be  used.  Hopefully,  the  Brady  sensor 
frequency  response  can  be  improved  or  corrected  for.  Part  of  the  frequency 
response  loss  may  be  due  to  the  protective  shield;  the  Brady  sensor  has 
been  relocated  to  improve  its  flushing. 

The  low  spectral  values  for  the  UBC  sensors  seem  peculiar  because 
visual  examination  of  the  original  records  shows  that  the  three  follow  one 
another  rather  well  for  periods  of  a minute  and  longer.  The  phase  difference 
of  about  20°  (Fig.  13b)  for  frequencies  about  0.01  Hz  would  suggest  that  the 
error  in  amplitude  would  be  moderate  if  an  RC  type  of  response  is  occurring. 
The  UBC  calibrations  have  been  checked  and  while  they  have  some  uncertainty 


the  large  discrepancy  cannot  be  accounted  for  by  calibration  errors;  the 
BIO  calibrations  must  be  reasonable  since  the  o /AQ  and  C values  are 

9 q 

reasonable . 

The  Brady  sensor  is  supposed  to  respond  to  relative  humidity.  In 


calculating  and  u ^q  values  temperature  effects  have  been  ignored  - i.e. 
the  relative  humidity  values  have  been  multiplied  by  the  interpolated 
saturated  values  of  Q based  on  the  air  temperatures  observed  with  the 
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TABLE  5 
HUMIDITY  DATA 


S = water  surface;  P = hand-held  psychrometer  at  9m; 

L = Lyman-alpha;  B = Brady;  A = aluminum  oxide 


Run 

Relative 

Humidity 

Water  Vapour  Density 

Q . AQ 

gm/nr 

Q 

& 

> 

qu3 

gm/m^s 

V1"3 

qu^/AQUio 

6 

S 

100 

12.1 

_ 

P 

68.7 

8.0 

4.1 

- 

_ 

_ 

L 

7.8 

4.4 

0.088 

0.050 

1.51 

B 

69.0 

8.1 

4.0 

0.026 

— 

— 

A 

7.5 

4.6 

o.oi3 

- 

- 

7 

s 

100 

12.1 

_ 

_ 

_ 

p 

69.7 

8.1 

4.0 

- 

- 

_ 

L 

8.1 

4.0 

0.097 

0.049 

1.66 

B 

10.0 

8.1 

4.0 

0.030 

0.007 

0.20 

A 

7.6 

4.5 

0.010 

0.006 

0.16 

8 

S 

100 

12.1 

_ 

_ 

_ 

_ 

p 

70.5 

8.3 

3.8 

- 

- 

— 

L 

8.3 

3.8 

0.091 

0.  047 

1.47 

B 

70.1* 

8.3 

3.8 

0.025 

0.004 

0.12 

A 

7.9 

4.2 

0.027 

0.006 

0.15 

9 

S 

100 

11.9 

_ 

_ 

_ 

_ 

P 

72.3 

8.8 

3.1 

- 

- 

_ 

L 

8.7 

3.2 

0.078 

0.026 

1.05 

B 

70.  4 

8.6 

3.3 

0.022 

0.002 

0.08 

A 

8.5 

3.4 

0.051 

0.003 

0.11 

C 

q 

= (1.42  i 

: 0.26)  x 

10~3 
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mercury  thermometers  of  the  psychrometer . The  temperature  fluctuations 
also  have  an  effect  but  under  the  observed  conditions  this  effect  is  only 
about  10$.  With  temperature  fluctuations  generally  in  phase  as  is  the 
case  for  these  runs  and  well  correlated  (as  is  usually  the  case)  the 
corrected  and  u^q  values  would  be  10$  larger  and  the  spectral  values 
20$  larger,  not  nearly  enough  to  account  for  the  observed  discrepancies. 

Both  the  Brady  and  AlO^  sensors  have  some  hysteresis.  Perhaps  the 
hysteresis  is  the  cause  of  the  low-frequency  attenuation  in  amplitude  with 
small  phase  shift.  The  higher  frequency  RC-type  rolloff  is  due  to  limited 
response  of  the  sensor. 

Water  surface  temperatures  were  measured  with  a bucket  thermometer 
and  the  saturated  air  in  contact  with  the  water  is  assumed  to  have  been  at 
the  same  temperature.  The  average  value  of  the  evaporation  coefficient, 

C = (1.42  ± 0.26)  x 10  for  the  Lyman-alpha  and  sonic  u data,  agrees  well 

3 

with  existing  observations  (e.g.,  Pond  et  at.,  1971;  Smith,  197*+;  Friehe  and 
Schmitt,  1976). 

5.  SUMMARY 

The  experiment  achieved  its  main  objective:  to  demonstrate  that 

the  BIO  system  and  the  IOUBC  system  would  run  well  together  and  measure 
essentially  the  same  quantities.  The  agreement  for  velocity  measurements 
was  good  if  corrections  were  applied  to  account  for  the  response  limitation 
of  the  Gill  anemometer  and  for  an  apparent  calibration  difference  between 
the  sonic  and  Gill  anemometers.  The  thrust  anemometer  results  were  more 
scattered  because  of  noisy  operation,  but  nevertheless  also  agreed  with  the 
sonic  and  Gill  anemometer  results  fairly  well.  The  overall  results  were 
consistent  with  similar  measurements  obtained  at  other  sites.  Temperature 
measurements  with  the  two  systems  showed  general  agreement  for  temperature 
fluctuation  and  fluxes,  but  not  for  mean  values.  There  were  problems  with 


ho 


calibration,  drift,  frequency  response,  and  noise  of  the  temperature  sensors 
which  we  shall  try  to  reduce  in  future  experiments.  The  limited  humidity 
measurements  showed  that  there  is  reasonable  coherence  between  sensors  at 
low  frequencies  but  that  the  IOUBC  sensors  have  reduced  response  at  moderate 
frequencies  as  well  as  a time  constant  effect  at  higher  frequencies. 

As  a result  of  this  trial,  development  of  a Mark  8 thrust  anemometer 
has  been  undertaken  for  measurements  at  a stable  platform  near  Halifax.  It 
will  use  proximity  transducers  instead  of  differential  transformers,  and 
requires  a less  elaborate  set  of  springs.  A counterbalanced  design  for  the 
horizontal  components  and  much  higher  resonant  frequency  ( > 60  Hz ) should 
substantially  reduce  the  vibration  problems  reported  here  for  the  Mark  6 and  7 
anemometers.  A thermistor  will  be  mounted  in  an  air  passage  in  the  case  of 
the  anemometer  with  its  socket  inverted  so  that  moisture  will  drain  out  instead 
of  in. 

The  need  for  a fast-response  humidity  sensor  which  can  survive  in  a 
marine  climate  remains  unfulfilled.  The  Lyman-alpha  sensor  performs  well 
wtien  carefully  handled  but  has  large  power  consumption  and  calibration  drift. 
The  Brady  sensor  needs  improved  frequency  response  and  calibration  stability. 
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